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We present quantum mechanical explanations for unresolved phenomena observed in field-free 
molecular alignment by a femtosecond laser pulse. Quantum phase analysis of molecular rotational 
states reveals the physical origin of the following phenomena: strong alignment peaks appear peri¬ 
odically, and the temporal shape of each alignment peak changes in an orderly fashion depending 
on molecular species; the strongest alignment is not achieved at the first peak; the transition be¬ 
tween aligned and anti-aligned states is very fast compared to the time scale of rotational dynamics. 

These features are understood in a unified way analogous to that describing a carrier-envelope- 
phase-stabilized mode-locked laser. 

PACS numbers: 37.10.Vz, 33.20.Sn, 33.80.-b, 42.50.Hz 


Chemical reaction depends on the relative orientation 
of reactants due to anisotropy in their electronic struc¬ 
ture. Likewise, in photochemical processes, the relative 
orientation between laser polarization and molecular axis 
is a basic parameter affecting the result. In order to clar¬ 
ify or to control the chemical processes at the most fun¬ 
damental level, the processes should be examined with 
aligned molecules. Among several techniques developed 
so far [ll-Q, the field-free alignment method using an in¬ 
tense femtosecond laser pulse is now widely adopted. Due 
to the anisotropy of molecular polarizability, the intense 
laser pulse induces a torque to rotate molecules to the 
direction of laser polarization. After interacting with the 
laser pulse, the molecules periodically align to that direc¬ 
tion, achieving field-free alignment [3, 0] . This method 
has become popular in revealing molecular structures 
controlling and tracing chemical reactions [3, Q, 
generating high-order harmonics from aligned molecules 
and compressing optical pulses [ll[ [T^ . 

The field-free alignment method exhibits unique fea¬ 
tures in the temporal variation of alignment. A typ¬ 
ical temporal structure of the method can be seen in 
Fig. 1, which is obtained by numerically solving the time- 
dependent Schrddinger equation (TDSE) for the rota¬ 
tional state of O 2 . The degree of alignment can be repre¬ 
sented by ((cos^0)), in which 9 is the angle between the 
laser polarization and the molecular axis and the double 
bracket means that the value is obtained by Boltzmann- 
averaging of all (cos^0) (quantum mechanical expecta¬ 
tion value) in a thermal ensemble. The unique features 
can be summarized as follows. First, it shows a revival 
structure with a full revival time Trev in which multi¬ 
ple fractional revivals occur, showing periodic alignment 
peaks Second, each alignment peak has its own 

temporal profile that changes in an orderly fashion de¬ 
pending on molecular species [l^. Third, the strongest 
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FIG. 1. (Color online) Temporal variation of the degree of 
alignment ((cos^^)) of O 2 molecules at 90 K after interacting 
with a sin^-type laser pulse with a duration of 30 fs (FWHM) 
and a cycle-averaged peak intensity of 5 x 10^® W/cm^. 

alignment is not achieved at the first peak [l^. Fourth, 
the transition between aligned and anti-aligned states is 
very fast compared to the time scale of rotational dy¬ 
namics [ 1 ^, as marked by ATtr in Fig. 1. These fea¬ 
tures have been observed also in experiments 0113 , and 
many theoretical studies have been made to understand 
the physical origin of the phenomena 0 Hill, BIS- 
Considering that the revival structure is a salient feature 
of the coherent superposition state of quantum systems, 
e.g. coherently excited Rydberg atoms [2l|, the complex 
amplitudes of the molecular eigenstates need to be ana¬ 
lyzed in detail to explain the features [Tsl - f^ . 

In this work we present an analysis of the molecular 
alignment induced by a femtosecond laser pulse in terms 
of the quantum phase of molecular rotational states. By 
numerically solving the TDSE of molecular rotational 
states, the quantum phase of rotational eigenstates is 
obtained. With the phase of the eigenstates known, the 
temporal variation of molecular alignment is described in 
a way similar to that of mode-locked laser pulses. The 
results provide the quantum mechanical explanation on 
all the above features in a unified way, signifying that the 
quantum phase of rotational states plays a crucial role in 
understanding the alignment dynamics. 

To be specific, the alignment of an O 2 molecule, a typ- 
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ical non-polar linear molecule, by a femtosecond laser 
pulse was considered for analysis. The TDSE for such a 
case within a rigid rotor model in atomic units is given 
by UM 

^ “ ^5(t)cos20^ , (1) 

where Bq is the molecular rotational constant, I the 
cycle-averaged peak intensity, Iq = Bqc/ (27rAQ;) the nat¬ 
ural intensity scale of the rotational molecular Hamil¬ 
tonian, and g{t) = sin^ {nt/r) the normalized laser- 
intensity profile. Aq;{= a\\ — a±) is the anisotropy of 
molecular polarizability. O 2 molecule has Bq = 2.856 x 
10“^® erg (1.437 cm“^) and Aa = 1.12 [ 2 ^, giving a 

natural rotational time scale tt/Hq = 11.6 ps and Iq = 
1.22 X 10^^ W/cm^. In our calculation, the condition of 
r = 60fs = 0.00527r/Ho and / = SxlQi^ W/cm^ = 4 IO /0 
was used. 

The evolution of molecular rotational states during the 
laser interaction is investigated by considering an angular 
momentum eigenstate \Ji,Mi) as an initial state. Since 
the laser interaction causes sequential Raman-type tran¬ 
sitions with the selection rules A J = 0, ±2 and AM = 0 
3 , the rotational state under interaction can be written 
as 


|4'(t))= ^ I + 2n, M), (2) 

n—no 

where Pj.+ 2 n and $j^_|_ 2 „ are the probability and the 
quantum phase of n-th emerging component, respec¬ 
tively. The emerging components with positive n come 
from Stokes processes, while those with negative n from 
anti-Stokes processes. Since Ji + 2n cannot be smaller 
than I Mi I, the anti-Stokes processes are bound by a 
non-positive integer uq [^ . Inserting Eq. (2) into 
Eq. (I) gives coupled differential equations of Pj^+ 2 „(t) 
and $j^_|_ 2 „(t). In obtain ing the numerical solution, the 
Crank-Nicolson method was used. Eigure 2 shows 
the numerical solution of the rotational state for an ini¬ 
tial state 13,0) as an example. 

The quantum phase $j^_|_ 2 „ is the basic element of 
the alignment dynamics. As seen in Fig. 2 (a), the 
quantum phase of n-th emerging component starts with 
TT |n|/2 [isl - l^ . This phenomenon can be understood 
by analyzing the mathematical form of Eq. (1). The 
emerging components are generated through the off- 
diagonal matrix elements of the interaction Hamiltonian, 
—Bolg{t)cos^9/ Iq, which are always negative as the ma¬ 
trix elements of cos^9 are positive [ 2 ^, and Aa is also 
positive for most linear molecules including O 2 . Because 
of this negativity of the interaction Hamiltonian and the 
factor i = exp(f7r/2) in the left-hand side of Eq. (1), the 
starting phase of an emerging component must add it/2 
to the phase of its parent component. Since a new com¬ 
ponent arises only through sequential transitions, 7r/2 is 
accumulated to its phase whenever a succeeding compo¬ 
nent emerges. Even when the populating process is of 
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FIG. 2. (Color online) Temporal evolution of the rotational 
state initially at | Ji = 3, Mi = 0) of O 2 during laser interac¬ 
tion: (a) phase and (b) probability. The laser condition is the 
same as that in Fig. 1. Time 0 and 60 fs correspond to the 
beginning and the end of the laser pulse, respectively. 


anti-Stokes (n < 0), the phase accumulation is also tt/2, 
not —TT 12, because the matrix element of the interaction 
Hamiltonian is equally negative. Thus, the n-th quantum 
phase starts at the value of tt |n| /2. 

Since the alignment is an interference among the ro¬ 
tational eigenstates, not the quantum phase itself but 
the relative phase determines the dynamics. As the laser 
interaction goes on, the quantum phase itself varies sig¬ 
nificantly, but the phase difference of 7r/2 between neigh¬ 
boring components set at the beginning does not change 
much because the phase-varying rates are comparable to 
each other, and the laser duration is very short compared 
to the time scale of the rotational dynamics. Thus, until 
the end of the laser pulse, the phase difference between 
neighboring components is kept to be around tt/2 [isl - l^ . 
as observed in Fig. 2 (a). 

After the interaction, the molecular wavefunction 
evolves freely and the dynamics of alignment is deter¬ 
mined solely by the molecular state at the end of the 
interaction (t = 60 fs in Fig. 2). Setting this moment as 
the new time origin, the field-free evolution of the degree 
of alignment, i.e. expectation value of cos^9, is given by 

00 

(cOS^e) (t) = ^ Pj,+2nCj,+2n,Ji+2u 

n—riQ 

+ |E (3) 

where 

= '^Pji+2nPji+2n+2ipJi+2n,Ji+2n+2) i ( 4 ) 

<^n ~ Pji+2n+2 ~ E'j._|_2ra = Bo{8n + AJi -|- 6), (5) 

— ^Ji+2n+2 ~ ^Ji+2n^ (6) 

Cj^j^ = {Ja, MI/\C0S^9 \ Jb, M/) , and Ej = BoJ{J + 1). 
Here, Pji+ 2 n and 4>j.+2„ are the values evaluated at the 
end of the laser pulse. Because Eq. (3) is expressed in 
a complex Fourier series, the time-domain behaviour can 
be analyzed in the frequency-domain, where the spectral 
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FIG. 3. (Color online) Spectral intensity (bars) and phase 
(filled circles) of (cos^6') (t) generated from several initial 
states with Mi = 0 after the laser interaction with molecules 
at 90 K: (a) Ji = 1, Mi = 0; (b) Ji = 3, Mi = 0; (c) Ji = 5, 
Mi = 0; and (d) Boltzmann-averaged spectrum from all the 
Ji states with Mi = 0 at 90 K. The laser condition is the 
same as that in Fig. 1. The illustrated spectral intensity 
was obtained by multiplying the corresponding Boltzmann 
weighting factor to Eq. (4). 


intensity frequency uin, and spectral phase (pn are 
given by Eqs. (4), (5), and (6), respectively. The first 
term of Eq. (3) is the DC term giving an offset in the 
time-domain, whereas the second term is the oscillating 
part having multiple frequencies. Equation (5) expresses 
an equally spaced frequency comb with a comb spacing 
of 8Bo. 

The spectral phase of (cos^0) (t) in Eq. (3) is the 
relative phase of two neighboring quantum states, as ex¬ 
pressed in Eq. (6). has nearly binary values close to 
either 7r/2 or —7r/2 [li,[1^, as shown in Eigs. 3 (a)-(c). 
This is explained by the initial quantum phases tt \n\ /2 
of which relative difference is mostly maintained until the 
end of the laser pulse. Because of this, Stokes processes 
(n > 0) result in ipn near tt/2 whereas anti-Stokes pro¬ 
cesses {n < 0) generate near —7r/2, making them nearly 
out of phase. This is the most critical part in under¬ 
standing the alignment dynamics. 

In order to obtain experimentally observable alignment 
dynamics, (cos^6>) (t) should be averaged over an initial 
thermal ensemble to give ((cos^0)) (t) as shown in Fig. 1. 
Instead of considering all the initial states at once, it is 
more systematic to consider only part of them by di¬ 
viding into groups according to Mi. Figure 3 illustrates 
a group with Mi = 0. For a given Mi, the lowest- 
state generates only Stokes processes so that the spec¬ 
tral phase is nearly constant, tt/2, as shown in Fig. 3 (a). 
Anti-Stokes processes produce low-frequency components 
and can happen only from initial high-states of which 
Boltzmann factors are smaller than that of the lowest- 
Ji state. Because of this, after cancelling out the spec¬ 
tral intensity of low-frequency parts due to the opposite 
phase between Stokes and anti-Stokes processes, the re¬ 
sulting ensemble-averaged phase follows the phase of the 
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FIG. 4. (Color online) Spectral intensity and phase of 
Boltzmann-averaged ((cos^S)) {t) from all initial states at 
90 K. The laser condition is the same as that in Fig. 1. 


Stokes processes as can be seen in Fig. 3 (d). In con¬ 
trast to the low-frequency parts, the spectral intensity of 
high-frequency parts add up constructively because only 
Stokes processes prevail in that region. Consequently, 
the ensemble-averaged phase of Fig. 3 (d) is formed al¬ 
most identical to that of Fig. 3 (a). Similarly, other 
groups having different Mi also show the same behav¬ 
ior, nearly constant phase of 7r/2. Owing to this phase 
coherence, the resulting spectral intensity from all the 
initial states is effectively enhanced. This is the main 
reason for achieving the strong alignment. The result¬ 
ing spectra of ((cos^6*)) [t) from all the initial states are 
shown in Fig. 4, where almost flat spectral phase and 
equally spaced frequency comb are observed. By syn¬ 
thesizing the spectral components in Fig. 4, the complex 
temporal variation in Fig. 1 was obtained. 

The complex temporal dynamics of molecular align¬ 
ment after the laser pulse can be understood by treating 
it as a mode-locked signal with almost flat phase. Since 
an equally spaced frequency comb having a linear spec- 
tral phase produces a pulse train in the time-domain | 26 l| , 
it is natural that ((cos^0)) (t) in Fig. 1 also shows pulsed 
shape like mode-locked laser pulses, exhibiting periodic 
strong alignment peaks. 

The temporal profile of each alignment peak in Fig. 1 
can be analyzed by using the analogy to the carrier- 
envelope phase (CEP) of a mode-locked laser pulse. In 
Fig. 1, the solid line corresponds to the carrier oscil¬ 
lation, while the dashed line the envelope. It is well 
known that the CEP-sequence is completely determined 
by the frequency-comb structure [^, [^. Let us de¬ 
note the frequency-comb as uim = ttiAlo + Slu {m = 
1, 2,---), where Aw is the pulse-repetition rate (mul¬ 
tiplied by 27r), and Suj is the offset frequency as shown in 
Fig. 4. Then, in the time-domain, the increment of CEP 
from one pulse to the next is Ac/cep = —Stt ■ Suj/Auj 
[^ . Since the full revival occurs when the accumulated 
CEP change reaches 27r, the number of fractional re¬ 
vivals occurring during a full revival period should be 
N = 2tt/ \A(I)cep\ = Aw/ |i5w|, and the full revival time 
should be Trev = 2ttN/ Alj = 27r/ |<5w|. The comb pa¬ 
rameters, Aw and 5lo, are obtained from Eq. (5) for 
a given molecular species. Since an O 2 molecule has 
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only odd-J states due to its nuclear spin statistics [^ . 
Aoj = 8 Bq and Suj = 2Bq as shown in Fig. 4, which 
results in A(j)cEP = — 7 '‘/ 2 , iV = 4, and Trev = t^IBq. 

The CEP of the first alignment peak is always nearly 
tt/2 (sine-like) regardless of molecular species, and it 
relates to the fact that the strongest alignment is not 
achieved at the first peak. The first CEP is determined 
by the temporal phases of all the frequency components 
at t = 0. These phases at t = 0 coincide with the spectral 
phases as can be deduced from Eq. (3). Since all the spec¬ 
tral phases lie near 7r/2, the first CEP should also be near 
7r/2. With the results of the last paragraph, the CEP 
sequence of O 2 molecule should be (7r/2, 0, —7r/2, —tt). 
This explains why the strongest alignment does not oc¬ 
cur at the first peak because it occurs when the CEP 
is 0. The same consideration gives the CEP sequence 
(7r/2, TT, — 7r/2, 0) for even-J molecules such as CO 2 m 
and {tt/2, —tt/2) for all-J molecules such as N 2 [P^ . 

Our analysis also provides a detailed explanation on 
temperature dependence of molecular alignment. It 
seems intuitively obvious that the degree of alignment 
decreases with temperature. However, although this phe¬ 
nomenon was partially explained from a classical treat¬ 
ment , no quantum mechanical explanation has been 
given. Here it can be explained by noting that Stokes and 
anti-Stokes processes generate opposite spectral phases. 
At a high temperature, considerable population is at 
high-Ji states initially. A higher-state causes more 
anti-Stokes processes that diminish the spectral intensity 
of low-frequency part. The higher the temperature is, 
the wider the frequency range is diminished. Because a 
small spectral intensity is equivalent to a low oscillating 
amplitude in the time-domain, the diminished spectral 
intensity results in weaker alignment. 

The fast transition between aligned and anti-aligned 
states can be explained by considering the center fre¬ 
quency of the frequency comb, ujc in Fig. 4. Wc cor¬ 
responds to the carrier oscillation period with the re¬ 
lation ATtr = tt/lOc- It means that high Wc is equiv¬ 
alent to fast transition. High cOc can result from ei¬ 
ther a high temperature or a strong laser interaction. 
A high temperature results in a high uic with two rea¬ 
sons. First, initially existing high-J^ states easily gener¬ 
ate high-frequency parts because, from Eq. (5), a;„ be¬ 
comes a high frequency with even a small n. Second, 
the spectral intensity of low-frequency part is diminished 


as explained before. A strong laser field also causes the 
high uic because the strong interaction populates much 
higher-J states through high-u Stokes transition so that 
the spectral intensity in the high-frequency region in¬ 
creases. Thus, ATtr becomes shorter when the tempera¬ 
ture is high or the laser interaction is strong. In Fig. 4, Wc 
is 29 . 8 H 0 leading to ATtr = 388 fs = 0.034Trev This fast 
transition is not surprising because the superposition of 
the eigenstates of J = 7 and J = 9 generates oJc = 34i?o 
leading to a transition time of 341 fs (= 0.030Ti.ev)- 

The molecular alignment at higher intensity was also 
examined. For example, at / = 8 x 10^^ W/cm^ = 655/o, 
the quantum phase $j. showed an abrupt jump due to 
the redistribution of rotational states during the laser 
interaction. Due to the phase jump, the tt/2 phase- 
difference between adjacent quantum phases did not hold 
any more, and the spectral phase of (cos^ 0 ) (t) deviated 
from that of Figs. 3 (a)-(c). It was, however, found that 
the Boltzmann-averaging mitigate these effects to make 
the ensemble-averaged spectral phase to be almost the 
same as that of Fig. 3 (d). 

In conclusion, we have presented clear quantum me¬ 
chanical explanations on the puzzling phenomena ob¬ 
served in the field-free molecular alignment. Analyzing 
the quantum phase of molecular rotational states, we 
have shown that ((cos^ 0 )) (t) has almost flat spectral 
phase and equally spaced frequency comb. All the fea¬ 
tures of complex alignment dynamics could be explained 
in an analogous way describing a CEP-stabilized mode- 
locked laser. The analysis given here can be extended 
further to more challenging subject such as orientation of 
polar molecules or three-dimensional alignment of asym¬ 
metric top molecules. The clear understanding for align¬ 
ment dynamics will bring the optimization of experimen¬ 
tal conditions and the development of a better method 
controlling molecular rotational states. 
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